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ABSTRACT 
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Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle 
Main Engine (SSME) and LOX flow distribution through the LOX posts have a great effect on 
the combustion phenomena inside the main combustion chamber. In order to design a CFD 
model to be an effective engineering analysis tool with good computational turn-around time 
(especially for 3-D flow problems) and still maintain good accuracy in describing the flow 
features, the concept of porosity was employed to describe the effects of blockage and drag force 

due to the presence of the LOX posts in the turbulent flow field around the main injector 

"f.>: ‘ , ' .■ 

assembly of the SSME. JMD ’numerical studies were conducted to identify the drag coefficients 
of the flows, both through tube banks and around the shielded posts, over a wide range of 
Reynolds numbers. Empirical, analytical expressions of the drag coefficient as a function of local 
flow Reynolds number were then deduced. The porosity model was applied to the turbulent flow 
around the main injector assembly of the SSME, and analyses were performed. The 3-D CFD 
analysis was divided into three parts" LOX dome, hot gas injector assembly, and hydrogen cavity. 
The numerical results indicate that the mixture ratio at the downstream of injector face was close 
to stoichiometric around baffle elements. 
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1.0 INTRODUCTION 


1.1 THE NATURE OF THE PROBLEM 

Fluid flow and heat transfer in the main injector assembly of the Space Shuttle Main 
Engine (SSME) are complex phenomena. The basic understanding these phenomena is essential 
to achieving optimum performance during normal operating conditions and maintaining structural 
integrity during off-design operations. The mixture ratio and mass flow rate distributions of the 
SSME main injector assembly will greatly affect 1 ) engine performance, and 2) heat loads of the 
combustion chamber; especially, the later effect is directly linked to the durability of the engine. 
Historically, the SSME has been suffering burn out of the LOX post baffle elements and erosion 
of the combustion chamber during firings. In order to investigate possible causes of such 
damage, understanding the flow field at main injector exit is essential. The geometry of the 
SSME main injector assembly is extremely complex, and its flow field is three-dimensional and 
turbulent. Conventional three-dimensional computational fluid dynamic (CFD) models are not 
practical to describe the necessary geometric detail of the SSME main injector assembly. The 
flow description was simplified by utilizing the concept of porosity to provide an effective 
engineering design tool for this system. 

1.2 OBJECTIVE 

The objective of this study was to develop a practical CFD simulation of the main injector 
assembly. The geometric complexity caused by the use of hundreds of individual LOX -post 
elements was reduced to a manageable computation by using non-isotropic porosity and 
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distributed resistance models. LOX-post arrays, shielded LOX-posts, and flows through porous 
injector plates were characterized with blockage and resistance models. A non-isotropic porosity 
model was incorporated into an existing Navier-Stokes flow solver (FDNS). Volume and surface 
porosity parameters, which are based on the configurations of local lox-post clustering, were 
introduced into the governing equations. Accuracy and robustness of the proposed model was 
demonstrated through data comparisons with benchmark test data and with detailed CFD 
solutions. Application of the postulated model to the turbulent flow within the main injector 
assembly of the SSME was made. This design tool predicts the local O/F distribution of the flow 
entering the main combustion chamber. 
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2.0 CFD METHODOLOGY 


The turbulent flow around the LOX post assembly is similar to the flow through a tube 
bank ensemble. There are basically three methods available in the literature (Ref. 1) to analyze 
the fluid dynamics and heat transfer characteristics of flow-cylinder assemblies with various 
configurations. These methods are: 1 ) sub-channel analysis; 2) porosity and distributed resistance 
approach (Ref. 2); and 3) benchmark detail rod-bundle fluid/thermal dynamics analysis using 
boundary fitted coordinate (BFC) system (Refs. 3,4). The first method is a simplified approach 
of the second method. Although the last method of analysis can provide the most detailed 
computational results, the mesh size required to resolve the geometrical complexity of the entire 
main injector assembly prohibits its use as an effective engineering design/analysis tool. One 
solution for this problem is the use of porosity modeling in the CFD analysis which will provide 
much better computer turn-around time. Validity of the approximations employed in the porosity 
model can be verified by comparing with the detailed CFD/BFC solutions for geometrically 
simplified test cases. A non-isotropic porosity model was developed and validated by comparing 
to the 2-D tube bank flows over a wide range of Reynolds numbers. The mass flow rate through 
the post element and the porous plate was also calculated by using the distributed resistance 
model. The non-isotropic porosity and the distributed resistance models were incorporated into 
the Navier-Stokes flow solver (FDNS). 
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2.1 GOVERNING EQUATIONS 

The FDNS code (Refs. 5-7) is a time-accurate pressure-based predictor-corrector flow 
solver. Various turbulence models, such as standard k-e mode, extended k-e model (Ref. 8), low 
Reynolds number k-e model, along with different compressibility corrections, have also been 
incorporated into the code. The FDNS code solves the following form of the conservation 
equations, including the continuity equation, the Navier-Stokes equations, the energy equation, 
and two-equation turbulence models, in curvilinear coordinates: 

1 3oa 3F. 

= - ■ +S q = R = Residual (1) 

J dt dfT q q 

where q stands for dependent variables (i.e. unity, the velocity vectors, temperature, turbulence 
quantities, and mass fractions of chemical species), and the numerical flux, F, is the sum of a 
convective flux, F c , and a viscous flux, F v , i.e. 


F = F + F 


Fc = 


pU.q 


F = -a G 

r v 


( 2 ) 


where J, U; and represent the Jacobian of the coordinate transformation, transformed 
velocities, and diffusion metrics, respectively. They are written as 


T _ dd.'H.O n - Ui r, - 1 n) 

d( x, y, z ) ’ 1 " J 3T ’ ij ■ Jd^d^ 

= (p, + p l )/o q is the effective viscosity when the turbulent eddy viscosity concept is 

employed to model the turbulent flows, p, = pC M k 2 /e, is the turbulence eddy viscosity and C M 

and a q denote turbulence modeling constants, and k and e are the turbulence kinetic energy and 
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its dissipation rate; while p, is the fluid viscosity. Source terms S, for the continuity, momentum, 
energy and species equations are given by 


0 
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where p and Uj are the static pressure and velocity vectors, V 2 = £ Uj 2 , P r is the turbulence kinetic 
energy production rate, co n is the species production rate, where C„ C 2 , and o T are turbulence 
modeling constants. The equation of state for an ideal gas or a real gas is used to close the 
above system of equations. These equations were solved with the pressure-based FDNS solver, 
which is discussed in detail in Refs. 5 & 6. 


2.2 NUMERICAL ALGORITHM 

A pressure based solution method is employed in FDNS so that a wide range of flow 
speeds can be analyzed with the same code. For high speed flow cases, a hyperbolic pressure 
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correction equation is retained by perturbing the density in the mass conservation equation. This 
provides a smooth transition from low speed to high speed flows. For time accuracy, a time- 
centered, time-marching procedure with a multiple pressure corrector algorithm is employed. 
This method provides numerical efficiency for time-dependent flow problems. 


To solve the system of nonlinear coupled partial differential equations, finite difference 
approximations are used to establish a system of linearized algebraic equations. A relaxation 
solution procedure is then employed to couple the equations. First, equation (2) is discretized 
in time with a time-centered (Crank-Nicholson) scheme. That is, 

^[(pqr'-ipqrl^^^i (5) 

where the superscripts n and n+1 denote the current time level and the next time level, 
respectively. If a sub-iteration procedure, designated by a superscript k, within a time step is 
applied, the following linearization is incorporated. 

( p q ) n *' = (pq) k + p n Aq k (6) 


. n+ 1 


f v 

3R, 

V 3c l j 


Aq k + R' 


(7) 


With the above linearization, the final form of the time-marching scheme is written as 


' p v 


f V 

dR 

q 


A< i l - - A-[(pq) t -(pq ) 11 ] * ■— 

J At 2 


( 8 ) 
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The solution at time level n+1 is then updated by 


q"* 1 = q k *‘ = q k + Aq k (9) 

when k = 1 is selected, a non-iterative time-marching scheme is used. As reported in Ref. 6, the 
non-iterative option with a multi-corrector solution method provides time accurate solutions for 
transient flow problems. This multi-corrector procedure is described briefly below. 

A simplified momentum equation is combined with the continuity equation to form a 
pressure correction equation. This pressure correction equation exhibits elliptical behavior for 
low speed flow and becomes continuously more hyperbolic as flow speed increases. The 
simplified momentum equation can be written as 

3pu a, 

-4-i * -Vp\ or u. « -^iVp' (10) 

3t 1 p 

where the superscript ' denotes perturbations. The velocity and density fields in the continuity 
equation are then perturbed to form a correction equation. That is 

V(p Ui r' = V[(p" + p / )(u i " + U i ')] (11) 

By neglecting the p'u/ terms, the following equation results. 

V(ap') + V(pu i / ) = -VCpUj)" (12) 

Equation (12) has the form of a transport equation with convection and diffusion terms. Upwind 
treatment can be used to model the left hand side of equation (12). A dissipation term is also 
added to the right hand side of equation (12) to provide smooth shock solutions. Using the 
solution of equation (12), the pressure and velocity fields are updated and the density field is then 
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updated by using the equation of state. This corrector procedure is repeated several time (usually 
4 times are sufficient) before marching to the next time step. This procedure insures that the 
mass conservation condition is satisfied for each time step. This represents the multi-corrector 
solution procedure. This method requires one predictor step and less than four corrector steps 
to provide numerical efficiency for transient flow computations. 


8 



3.0 PROPOSED POROSITY MODEL 


The conventional porosity models assume the use of orthogonal coordinates and 
geometrically similar control volumes (Refs. 1 & 2). In the proposed approach, the general 
boundary fitted coordinate systems were incorporated in the formulation. Two new parameters, 
volume porosity (y v , defined as the ratio of the volume occupied by the fluid to the total volume) 
and surface porosity (y, defined as the ratio of the surface area in the i-direction available for 
flow to the corresponding total surface area in the same direction), were introduced into the flow 
governing equations. For the proposed porosity model, equation (1) was rewritten as 


y v 9pq 

TIT 




+ y S 

• v q 



= Residual 


(13) 


A distributed drag force D and a heat flux H source term were also added to the right 
hand side of the momentum equation and the energy equation, respectively, to simulate the 
effects of resistance and heat transfer due to the presence of the LOX posts in the flow field. 
These drag force and heat flux terms were modeled based on geometric parameters and the 
averaged velocity around a local LOX post. Since the drag force D is defined as 


D = 0.5p u 2 C d (14) 

where p, u, and C D are local flow density, local total velocity, and local drag coefficient, 
respectively, we can compute the distributed drag force by evaluating these three parameters. 
2-D tube bank flow studies were investigated to verify the drag coefficients for the flow around 
the LOX posts assembly. Moreover, the 13th row posts are shielded in pairs, where on the shield 
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surface there are several holes distributed between each pair of posts. Additional CFD validation 
studies were needed to identify the drag coefficients for the flow around shielded elements. 

3.1 2-D TUBE BANK FLOW STUDIES 

The CFD investigation of the flow through the tube bank configuration was conducted 
for various Reynolds numbers (Re), such as 20, 10\ 10 6 , and 10 7 . The flow fields of the Re=20 
case and the Re=10 6 case are demonstrated in Figures 1 & 2, and the computed pressure 
distribution along the cylinders with these four Reynolds numbers are plotted as shown in Figures 
3-6. 
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Figure 1 Streamline Velocity (u) contours for the Flow Past a 2-D Tube Bank with Re=20 
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Figure 2 Streamline Velocity (u) contours for the Flow Past a 2-D Tube Bank with Re=10 6 


10 





1 0E*02 


9 9E*01 

P/C2Q) 

9. 6E + 01 


9 3£+0 1 


9 0E+0 1 1 

• M 


Figure 5 


Static Pre 


i 0E+02 r 


9 9E + 01 


P/C2Q) 


9.6E+01 


9 3E+01 


9 0E+01 

IK 


Figure 6 Static 





Based on the CFD studies and a literature review (Refs. 9-14), a drag coefficient model 
for the tube bank flows is proposed as shown in Table 1. The proposed model was compared 
with the numerical results from previous benchmark cases and relatively good agreement has 
been achieved. 


Table 1. The Drag Coefficient Model for Tube Bank Flow 


1 

C D 

•!'>Xv!'X*X'fi*X*XvIv!v'X\5/vS® 

Re < 4 x 10 3 

0.417 EXP(4.932 Re- 0 - 296 ) 

4 x 10 3 < Re < 6 x 10 4 

0.647 - 0.5 x 10’ 6 Re 

6 x 10 < Re < 10 6 

0.618 + 0.491 x lO^Re - 6.303 x 10 12 Re 2 
+ 10.694 x 10 l8 Re 3 - 5.2 x 10 24 Re 4 

Re > 10 6 

0.2735 


3.2 TURBULENT FLOW AROUND SHIELDED POST ELEMENTS 

In the hot gas injector region, the outer row (#13) is protected by shields to avoid 
damages caused by direct impingement on the LOX posts from the high speed gases. There are 
four types of shield, #039, #037, #025, and #023, where the configurations of shield-039 and 
shield-023 are sketched as shown in Figure 7. All the shields enclose two posts in one shield, 
except shield-025 which encloses three posts in one shield. The hole distribution on shield-025 
is similar to that on shield-023, and there is only one shield-025 on row #13; hence shield-025 


13 


Computational Domain 



Figure 7 The Configuration of the Various Shielded Elements 
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was modeled the same as shield-023 in the numerical analysis. Shield-037 has the same type of 
hole distribution as shield-039, but shield-037 does not have the wing at the end of the shield. 
There are six #037 shields and five #039 shields, and they are installed in an alternating sequence 
and are located near the fuel side. In the proposed porosity, the drag coefficient for both shield- 
039 and shield-037 was assumed to be identical and was obtained from the numerical analysis 
of shield-037 geometry. Therefore, in this study all shields were considered to be of either the 
023 or 037 type. CFD analyses were conducted for these two types of shields to obtain the 
porosity model for the flow through a shielded element. 

Since the presence of the shield will greatly increase the drag force, numerical studies 
were performed to investigate the drag coefficient for the flow through the shielded elements with 
and without holes. Both 3-D and 2-D analyses were performed. Since the shielded elements are 
located on the outer row of the post assembly, the Reynolds number of the inlet flow to the 
shielded elements is relatively high and is around 10 6 . A symmetry boundary condition was 
specified at the centerline of the gap between posts to minimize the computational domain as 
indicated in Figure 7. The effect of the T-bolt as shown in Figure 7 was not taken into account 
in this study. An approximation for different types of holes on the shields was made to relax the 
requirement of using numerous grids to resolve these holes. The approximation was to treat the 
cross-section of various holes as rectangles with their width equal to the gap between each pair 
of posts. Therefore, the same drag coefficient model would be employed for various types of 
shields. The way to distinguish the difference between various types of shield is to identify the 
vertical locations of each hole and to use different porosity for different types of shields. 
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The numerical studies indicate that the drag coefficient does not change appreciably as 
the Reynolds number varies around 10 6 . The computed drag coefficient is close to 4 for the 
region with holes on the shield, and is close to 48 for the region without holes. The calculated 
flow fields for the flow through the shielded element with holes and without holes are exhibited 
in Figures 8a & 8b, respectively. 



Figure 8a Pressure Contours and Velocity Vectors for the Flow Past a Shielded Element with 
Holes 



Figure 8b Pressure Contours and Velocity Vectors for the Flow Past a Shielded Element 
Without Holes 
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3.3 LOSS COEFFICIENT MODEL 

In addition to the drag coefficient model for the flow through tube bank environment, a 
loss coefficient model is required to compute the mass flow rate for the flow through various 
porous media, such as an individual LOX post, porous plates, and boundary layer control (BLC) 
holes. To relate the mass flow rate and pressure drop across the porous media, the following loss 
coefficients ( K) are defined. 

K = , Ap = p - p ch . mber , or Ap = p - Pbifne (15) 

nr 

where p, p are the density and static pressure at the exit plane, and m is the mass flow rate 
through a porous medium. The value of the loss coefficient ( K ) for each porous medium will 
be given later in the NUMERICAL RESULTS section. Therefore, the exit mass flow rate 
through each porous medium was computed based on the distributions of local loss coefficient 
and local exit pressure. In the numerical analysis, the area for each porous medium is different 
at each radial location; however, this is not true in the real geometry. Hence, in the 
computational domain, the area effect was excluded from the mass flow rate calculation. The 
assembly of the loss coefficients in the computational domain is then defined as 

(16) 

where K 1 and M are the global loss coefficient and mass flow rate for a type of porous medium, 
respectively. Hence, the local mass flow rate m can be calculated as 
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. A. 

m = M _L = 

A N 


P Ap A , 
K' T 


( 17 ) 


where A and A; denote the global and local exit areas for the same type of porous medium, 
respectively. The value of the global loss coefficient ( K *) for each porous medium will also be 
reported later in the numerical results. 
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4.0 NUMERICAL RESULTS 


The 3-D porosity/CFD analysis of the SSME main injector assembly was performed for 
the phase 11+ power head geometry at the 104% power balance level. The entire configuration 
was divided into three components: 1) LOX dome, 2) hot gas injector assembly, and 3) hydrogen 
cavity. The overall mass flow rate distribution and mixture ratio distribution was calculated by 
superimposing the results from these three components. The numerical computation for these 
three components was conducted independently, except the calculated exit pressure of the hot gas 
injector assembly was employed as the back pressure for the secondary face plate in the hydrogen 
cavity region. The exit pressure to the combustion chamber p chambcr was assumed to be uniform 
and equal to 3135 psi. The exit pressure from baffle elements p bafflc was also assumed to be 
constant and have the value of 3084 psi. 

4.1 LOX DOME 

The LOX dome geometry was simulated as shown in Figure 9 with a 62 x 91 x 16 mesh 
system. The bleed pipes to the first three rows of the LOX posts were omitted and were replaced 
with the exit flow through the nose region. This simplification reduced the number of grid points 
required and required less detailed information for the bleed pipe geometry. The importance of 
the mass flow rate distribution for the first three rows of the LOX post was mitigated by the 1 -D 
analytical result which indicates that the magnitude of the mass flow rate is relatively uniform 
among these three rows. The numerical analysis was conducted based on incompressible, 
adiabatic, turbulent flow with single species (liquid oxygen). The inlet flow conditions and loss 
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Figure 9 The Mesh System for the LOX Dome Geometry 






coefficients are listed in Table 2 and 3, respectively. Uniform inlet velocity profiles were used 
because at the LOX tee exit no other information was available. 


Table 2 Inlet Flow Conditions to the LOX Dome (104% RPL) 


Static Pressure 

Static Temperature 

Reynolds Number 

Mass Flow Rate 

3670 psi 

197 °R 

1.28 x 10 8 ft' 1 

826.7 lb/sec 


Table 3 The Loss coefficient model for the LOX dome region (104% RPL) 



m (lb/sec) 



HII8 1 

Non-Baffle Elements 

1.38 

9.62 x 10 7 

665.105 

4.14 x 10 2 

Baffle Elements 

1.433 

9.78 x 10 7 

105.65 

1.79 x 10 4 

First Three Rows 

1.332 

1.03 x 10 8 

56.154 

5.81 x 10 4 


Based on the above inlet flow conditions and the proposed loss coefficient model, a 3-D 
numerical computation was performed. Figure 10, top view of the velocity vector plot, exhibits 
the LOX flow around the LOX dome. The cross-sectional views of the four velocity vector 
planes in Figure 10 are plotted as shown in Figures 1 la-1 Id. It is clearly shown that two 
recirculation zones occur in the plane with inlet flow, but a very smooth flow structure is formed 
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Figure 11a Velocity Vectors in the LOX Dome at the Cross Section of -90° (Fuel Side) 



Figure lib Velocity Vectors in the LOX Dome at the Cross Section of -2' 
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Figure 1 lc Velocity Vectors in the LOX Dome at the Cross Section of 90° (Oxidizer Side) 



Figure lid Velocity Vectors in the LOX Dome at the Cross Section of 178' 
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in the plane without inlet flow. It can be seen that the exit flow velocities are fairly uniform in 
each cross section. The static pressure contours at the LOX dome exit plane are drawn in Figure 
12. Although it seems higher pressure around the planes with inlet flow, the overall pressure 
difference is very small, which indicates a fairly uniform exit pressure. The contours of the exit 
velocity from the LOX dome to the injector face are plotted as shown in Figure 13, which 
demonstrate the same characteristics as that of the exit pressure. 

4.2 HOT GAS INJECTOR ASSEMBLY 

The computational domain for the hot gas injector assembly consists of three zones: 1 ) 
LOX post assembly torus, 2) fuel transfer duct, and 3) oxidizer transfer duct. A 37 x 91 x 25 
mesh system was employed for the LOX post assembly torus, and a 10 x 21 x 17 mesh system 
was used for the fuel transfer duct, while the oxidizer transfer duct was described by a 10 x 15 
x 15 grid system. The geometry and the grid system of the hot gas injector assembly is plotted 
as shown in Figure 1 4. The inlet flow conditions to both transfer ducts are tabulated in Table 
4. The inlet velocity and pressure profiles of the hot gas flow to the fuel transfer duct were 
interpolated based on the numerical results by Yang, and his coworkers (Ref. 15), while the inlet 
profiles to the oxidizer transfer duct were interpolated based on the image of Yang's results. The 
secondary velocity vectors of the inlet flows to the fuel transfer duct and the oxidizer transfer 
duct are sketched as shown in Figures 15a & 15b. 

The loss coefficient distribution used in the LOX post assembly was based on the air flow 
test data for various post elements. The value of the loss coefficients is listed in Table 5. 
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Figure 12 The Contours of Exit Pressure in the LOX Dome (psf) 
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Figure 13 The Contours of Exit Velocity in the LOX Dome (ft/sec) 


Oxidizer Side 
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Figure 14 The Geometry and Grid System of the Hot Gas Injector Assembly 



Figure 15a The Secondary Velocity Vectors of the Inlet Flow to the Fuel Transfer Duct 



Figure 15b The Secondary Velocity Vectors of the Inlet Flow to the Oxidizer Transfer Duct 
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K is the measured loss coefficient for each single post, while K' is the loss coefficient for each 
type of the posts in the computational domain. According to Equations (15-16), the relation 
between K and K! is 


K , _ pAp _ pAp _ K_ 
M 2 (Nm) 2 N 2 

where N is the number of each type of post elements. 


Table 4 The Inlet Flow Conditions to the hot gas injector Assembly (104% RPL) 


1 

Nr 


(fr*) 

Rate (lb/sec) 

O/F Ratio 







Fuel Side 

3351 

1666 

3.17 x 10 7 

77.55 (in half) 

0.8685 

LOX Side 

3353 

1254 

7.88 x 10 7 

33.38 (in half) 

0.599 


Table 5 The Loss Coefficients for the Hot Gas Flow Through the LOX Post 



Elements 

Baffle- 

Elements 

Row #13 

Row #12 

Row #1 - #11 

K (in' 4 ) 

135 

156 

152 

oo 

K* (ft 4 ) 

4.374 x 10 2 

5.751 x 10 2 

23 

oo 
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The numerical analysis for the hot gas injector assembly was conducted based on 
symmetrical, incompressible, isothermal, turbulent, no-reacting flow with two species (hydrogen 
and oxygen). Special attention was paid to the distribution of the porosity and of the drag 
coefficient on the 13th-row elements to identify the presence of holes on the shield. The 
predicted velocity vectors at the symmetrical plane are drawn as shown in Figure 16. It appears 
that the hot gas is deflected by the shielded elements (the outermost row), and passes through 
the non-shielded area of the posts. In addition, the exit velocity of the hot gas seems to be very 
uniform, except at the baffle elements through which there is no hot gas exit to the injector face. 
The deflected hot gas not only passes through the non-shielded area of the posts, but also sweeps 
around the torus, which can be seen from the velocity vector plots as shown in Figures 17-18. 
The hot gases from the fuel and the oxidizer sides flow into the region between the LOX posts 
and out into the main combustion chamber. The fuel side gases penetrate much further into a 
horizontal plane through the LOX-post region. Significant mixing occurs where the hot gases 
from the two sides of the engine meet. The exit velocity contours of the hot gas through the 
injector face are plotted in Figure 19. This figure shows that the exit velocity of the hot gas is 
somewhat larger near the fuel side than near the oxidizer side; however, the difference is very 
small which indicates that the flow is nearly uniform except at the baffle elements. The hot gas 
exit pressure contours on the bottom surface of hot gas injector assembly are plotted as shown 
in Figure 20, which exhibit the same characteristics as that of the exit velocity. The mixture ratio 
of the hot gas exit to the injector face is uniformly increased from the fuel side towards the 
oxidizer side as indicated in Figure 21, where the discontinuity at the symmetry plane is due to 
the linear extrapolation from interior points. 
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Figure 16 The Hot Gas Row Velocity Vectors at the Symmetry Plane of LOX Post Assembly 
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Figure 17 The Hot Gas Flow Ve 
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Figure 18 The Hot Gas Flow Vel 




CON FOUR IFUFLS 



aoaoaooo 


aaaacan : =j a 


i a a a a a a 


gQoaagoaoaoQOQQac; -aaaaaaaaaa 

oaoaaaaaaaooaoac- r. naaaaaaaaa 

SCDOlQ'-NnTinCONXtJia’-N:-' 

NNNnnnnncinnnnTTT" -^TTT^unnjiuij^ 


36 


ORIGINAL PAGE 
COLOR PHOTOGRAPH 


Figure 20 Hot Gas Exit Pressure Contours of the LOX Post Assembly (psf) 
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Figure 21 Hot Gas Exit O/F Ratio Contours of the LOX Post Assembly 


4.3 HYDROGEN CAVITY 


The hydrogen cavity region was descritized into a 29 x 91 x 14 mesh system, and is 
sketched as shown in Figure 22. The inlet flow conditions are listed in Table 6 based on 104% 
power balance. The inlet velocity profile was assumed to be uniform everywhere, and all 
velocity vectors were directed towards the center of radius. 


Table 6 Inlet Conditions of the Fuel into Hydrogen Cavity Region (104% RPL) 


Pressure (psi) 

Temperature CR) 

Reynolds no. (ft*).f§ 


3395 

449 

2.52 x 10 7 

14.55 (in half) 


The loss coefficient model is listed in Table 7, where the mass flow rate through each 
porous medium was based on 1-D analytical results. The back pressure of the secondary face 
plate was set to be 3288 psi from the hot gas injector assembly numerical result. An assumption 
was made that all hydrogen flow exiting through the secondary face plate to the hot gas region 
passed through the non-baffle elements and exited at the injector face. Therefore, in the post 
processing, the mass flow rate through the secondary face plate was added to the mass flow rate 
through the primary face plate at the same grid location. 

The 3-D numerical analysis for the hydrogen cavity region was performed based on 
symmetrical, incompressible, adiabatic, turbulent flow with single species (hydrogen). The 


38 





O 


39 


Figure 22 The Mesh System for the Hydrogen Cavity (29 x 91 x 14) 



numerical result, as shown in Figure 23, indicates that a recirculation zone occurs as the fuel flow 
enters the hydrogen cavity through a narrow gap. It also reveals that the velocity of the flow exit 
through the baffle elements is quite uniform, and is much larger than that exit through porous 
plates. The flow exit through BLC holes has the largest velocity, but the exit mass flow rate is 
relatively small due to small exit area. The uniformity of the exit fuel flow can also be seen 
from the contours of the exit flow velocity and of the exit pressure in Figures 24 & 25, 
respectively. 


Table 7 The Loss Coefficient Model Used in the Hydrogen Cavity Region (104% RPL) 



m (Ib/sec) 



Primary Face Plate 

6.77 

251 

3.578 x 10 4 

Secondary Face Plate 

3.41 

98 

5.506 x 10 4 

Baffle Elements 

15.25 

301 

8.467 x 10 3 

Non-Baffle Elements 

0 

251 

oo 

BLC Holes 

3.67 

251 

1.16 x 10 5 
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Figure 23 The Velocity Vectors of the Fuel Flow at the Symmetry Plane of the Hydrogen Cavity 


Hydrogen Caoity, Exit FIom Uelocity Contours (ft/sec) 
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Figure 24 The Contours of the Exit Flow Velocity in the Hydrogen Cavity (ft/sec) 
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Figure 25 The Contours of the Exit Pressure in the Hydrogen Cavity (psf) 


4.4 MAIN INJECTOR ASSEMBLY 


After conducting the numerical computation for the three components of the main injector 
assembly, a post processing calculation was made in order to obtain the overall O/F ratio at the 
injector face. By superimposing the exit mass flow rates and the mass fractions of hydrogen and 
oxygen from the three components of the main injector assembly, the mixture ratio distribution 
at the injector face was obtained as shown in Figure 26. It appears that the highest O/F ratio 
occurs near the baffle elements. The O/F ratio distribution also exhibits higher mixture ratio near 
the fuel side than that near the oxidizer side. The O/F distribution in the circumferential 
direction is of great interest, especially near the outer edge of the injector face, and so it is 
plotted as shown in Figures 27 & 28, respectively. Figure 27 plots the value without adding BLC 
coolant flow in the O/F ratio calculation, while Figure 28 shows the O/F ratio with BLC coolant 
flow added. Each spike in the mixture ratio plots occurs at baffle elements. The O/F ratio 
distribution in the radial direction is displayed in Figures 29 & 30, which shows the mixture ratio 
at the plane of -90° (on the fuel side) and 90° (on the oxidizer side), respectively. The plots 
reveal that the O/F ratio is close to stoichiometric around baffle elements. 
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Figure 26 The O/F Ratio Contours of the Main Injector Assembly at 104% Power Level 
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Figure 27 
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Figure 28 
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Figure 29 
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Figure 30 The 



5.0 CONCLUSIONS AND RECOMMENDATIONS 


The 3-D CFD/porosity analysis reveals that the mass flow rate at the injector face is 
relatively uniform. The predicted mixture ratio is close to stoichiometric (O/F ■= 8) around baffle 
element at 104% power balance level, which might cause very hot spots around baffle elements. 
However, due to the many assumptions made in the porosity model and the use of a very coarse 
grid system, the numerical results can only provide a qualitative trend. As can be seen from the 
loss coefficient model, the local mass flow rate distribution is dependent on both the pressure 
drop and the loss coefficient distributions. Hence, the assembly of the loss coefficient model is 
critical to the numerical result, and the availability of the measured loss coefficients for each 
porous medium will greatly improve the CFD analysis. In addition, the distribution of chamber 
pressure and of baffle element discharge pressure was assumed to be uniform in this study; 
hence, the CFD/porosity model can be improved if the actual distribution of discharge pressure 
were known and specified in the calculation. Meanwhile, a proper inlet flow profile to the LOX 
dome and to the hydrogen cavity can help the developed model to predict the flow field more 
accurately. The developed CFD/porosity model should be further tested at different power 
balance levels. The numerical results of this study should be used as the inlet conditions to the 
combustion chamber in order to predict the engine performance and heat loads of the chamber. 
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